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Oh; 

' We discuss various tests of the factorization hypothesis making use of the close 

relationship between semi-leptonic and factorized nonleptonic decay amplitudes. It is 
. pointed out that factorization leads to truely model-independent predictions for the 

ratio of nonleptonic to semi-leptonic decay rates, if in the nonleptonic decay a spin 
^ ' one meson of arbitrary mass or a pion take the place of the lepton pair. Where the de- 

' cay constants of those mesons are known, these predictions represent ideal tests of the 

. factorization hypothesis. In other cases they may be used to extract the decay con- 

G^ ! stants. Currently available data on the decays D^tt^ , D*+7r^, , D*^ 

I are shown to be in excellent agreement with the factorization results. A weighted 

^ I average of the four independent values for the QCD coefficient oi extracted from the 

^ ' data gives ai = 1.15 ± 0.06 suggesting that it may be equal to the Wilson coefficient 

ci(/i) evaluated at the scale = rrib. 

Oh! 
> 

X 

' The dynamics of nonleptonic weak decays is strongly influenced by the confining color 

forces among the quarks. In contrast to semi-leptonic transitions, where the lepton current 
naturally factorizes and one is left with the hadronic matrix element of a color-singlet quark 
current, nonleptonic processes are complicated by the phenomenon of quark rearrangement 
due to the exchange of soft and hard gluons. The theoretical description involves matrix 
elements of local four-quark operators, which are much harder to deal with than current 
operators. 

A great simplification can be accomplished if one is willing to adopt the factorization 
hypothesis, which relates the complicated nonleptonic decay amplitudes to products of me- 
son decay constants and hadronic matrix elements of current operators similar to the ones 
encountered in semi-leptonic decays. Despite its remarkable success in the description of 2- 
body decays of B- and D-mesons, precise tests of the factorization hypothesis are of utmost 
importance in order to find out its realm of applicability as well as its limitations. While 
many tests have been suggested or already carried out most of them do not simply 

test the factorization hypothesis, but rather factorization together with some phenomeno- 



1 



logical model or, alternatively, together with heavy-quark symmetry for dealing with the 
hadronic current matrix elements. It is the main objective of this short note to concentrate 
on such tests that do not suffer from additional uncertainties due to our unsatisfactory ways 
of dealing with non-perturbative QCD. 

As there exist several versions of factorization in the literature, let us begin by giving 
an unambiguous prescription of how to calculate the rate of some exclusive nonleptonic B- 
decay in the factorization approximation. We will concentrate on 6 — > c transitions, which 
are induced by the effective Hamiltonian 



G 



F 



Heff — —7= Vcb 



Cl(/i) Qt + C2ifi)Qf 



V2 

+ penguin operators. (1) 

It consists of products of local four-quark operators with scale-dependent Wilson coefficients 
Cj(/i). The operators Qi and Q2, written as products of color-singlet currents, are given by 

Q2 = {cu),-A {d'h),-A + {cc),-A {s'h),-A , (2) 

where d' and s' denote weak eigenstates of the down and strange quarks, respectively, and 
{ch)^_^ = c7/i (1 — 75) & etc. The Wilson coefficients of so-called penguin operators |^ 
in Eq. (P are very small. Their contribution to the dominant decay amplitudes may be 
neglected. 

If QCD was turned off, the Wilson coefficients of the operators Qf' and would be 
Ci = 1 and C2 = 0. These values are modified by hard gluon exchange. Evaluated at the 



scale fx = rrih ^ 5.0 GeV one finds in leading logarithmic approximation |T0[: ci(mb) = 1.12 
and C2{mb) = —0.26. 

According to the factorization hypothesis one may now write the hadronic matrix el- 
ements of Qf and as products of two current matrix elements As an example. 



we consider the decay amplitude of the transition n , which in the factorization 

approximation is given by 

^fact ^ ^ y^^ y.^ ^ ^- | ^ | ) ( | (c 6)y | B° ) . (3) 

Class I transitions like the one considered above, in which only a charged meson can be 
generated directly from a current, are proportional to the QCD coefficient ai. Its relation 
to the Wilson coefficients will be discussed below. Correspondingly, those decays in which 
the meson generated directly from the current is neutral, like the J/\&-particle in the decay 
B K J/^, are called class II, and their decay amplitudes are proportional to the QCD- 
coefficient 02. Factorized amplitudes in which there is interference between ai- and a2-terms 
are categorized as class III. 

Usually, form factor suppressed weak annihilation topologies (IV-exchange and quark- 
annihilation diagrams) are neglected in the calculation of factorization amplitudes. This is 
not an inherent property of the factorization approximation. Rather it is necessary from a 
practical point of view, since little more is known about form factors at such large time-like 
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momentum transfer than that they should be strongly suppressed. What really is an inherent 
property of the factorization approximation is the neglect of final state interactions (FSI). 
However, unlike Z)-decays, the decays of 5-mesons do not take place in a resonance region. 
Thus one has good reason to believe that ignoring the effects of FSI is a good approximation 
in 5-decays. 

Let us now turn to the relation between the Wilson coefficients and the QCD coefficients 
ai and 02. Naively, one would expect ai = ci(/i/) +^C2{f^f) and 02 = C2(/i/) +.^ci(/i/), with 
^ = 1/Nc, and fif denoting the factorization point, in 5-decays usually identified with mi,. 
However, experience in D-decays has shown that setting ^ = allows for a better description 
of the data, and it has been suggested to treat ^ or even ai and 02 independently as a free 
parameters |jl2|. Thus one can test the factorization hypothesis by checking whether or 
not the values for the QCD coefficient ai (02) as extracted from different class I (class II) 
transitions agree with each other. For ai also an absolute prediction becomes possible, by 
observing that varying the parameter ^ in the range < ^ < 1/3 induces no more than a 
10% change in ai. One would therefore expect ai = 1.1 ± 0.1 which has been confirmed 
in a recent extraction of ai from all available nonleptonic 2-body decays of 5- mesons [0. 
However, it should be stressed that the fit for ai has been strongly dominated by the two 
decay modesQ 5° ^ tt". It remains to be seen whether the same value of ai will be 

found from other decay modes as more precise data become available. 

As first pointed out by Bjorken, the close relationship between factorized amplitudes and 
semi-leptonic decay amplitudes provides the most direct test of the factorization assump- 
tion To this end, a nonleptonic decay width is related to the corresponding differential 
semi-leptonic decay width evaluated at the same g^. Let us consider the ratios 
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where fp is the decay constant of the pseudoscalar meson P, Vij is the appropriate KM- 
matrix element (associated with P) and (in the limit of vanishing lepton mass) 
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The helicity amplitudes Hoi^q^) and H±{q^) 
Bjorken has suggested this test with P 
less than 0. 

For heavier pseudoscalar mesons, Xp and X*p become model dependent and may quite 
substantially deviate from 1. In the infinite quark mass limit, one finds, for example, Xp^ ~ 



TT, in which case ~ X* ~ 1 to within 
as can be shown by expanding those quantities in powers of m^/ml 0. 



1.36 and X 



Ds 



0.37 



^Here and in the following, "!)(*)" stands for "Z? or D* 
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We can get rid of this model dependence altogether by replacing the pseudoscalar meson 
P in Eq. (|^) by a vector or pseudovector meson. In the factorization approximation one 
then finds 

4 = , r^.^„ ^ . = 67rVy lair \V^^' , (6) 



where now, of course, the KM-matrix element is associated with the (pseudo-) vector meson 
and fv denotes its decay constant. The reason for all form factors and additional kinematical 
factors to cancel in the ratio can be easily understood. For zero lepton masses, the lepton 
pair that in the semi-leptonic decay is generated by the {V — A) current carries spin one 
in its cm. frame. Integrated over the lepton angles keeping = [pi + p,^)'^ fixed, the 
production of the lepton pair is therefore kinematically equivalent to the production of 
a (pseudo-) vector particle with four-momentum (summed over all polarizations of the 
(pseudo-) vector particle). Corrections to Eq. (H) due to finite lepton masses are of order 
mj/my. With the ^-meson being the lightest spin-one meson these corrections may safely 
be neglected for electrons and muons. 

Setting V = g, we can use Eq. (^) to obtain two independent values for ai, since the 
decay constant fg is known^. These values should be compared with those obtained from 
Eq. (|^) with P = TT. However, as long as the differential spectrum of the semi-leptonic 
decay £^ v has not yet been measured, we must again resort to some form factor 

model in decays with a D-meson in the final state. In Table |, we have used the predictions of 
Ref. 0] for those two decays. They are based on an Isgur-Wise function extracted from data 
on the decay D*^ v with perturbative QCD-corrections and (model-dependent) 

1 / mg-corrections added on. Nonleptonic decay data used in Table |I| has been taken from 
CLEO fl^ and ARGUS [|15|. The ARGUS data as well as the predictions of Ref. have 



been rescaled using the new CLEO measurement BR(D*+ — > D'^vr^) = (68 ± 2)% [|T6 
The experimental number for the ratio where some of the systematic uncertainties 

drop out, has been taken from CLEO alone. We observe that all four values for the QCD 
coefficient a\ presented in Table | are in excellent agreement with each other and with the 
expectation from perturbative QCD, thus providing strong support for the factorization 
hypothesis in i?-decays with large recoil. Taking the weighted average of all four values 
gives a\ = 1.15 ± 0.06, which suggests that, just like in D-decays, we may have ^ = 0, i.e. 
the QCD coefficient ai may be equal to the Wilson coefficient Ci(/i) evaluated at the scale 
of the decaying quark. 

As better statistics becomes available the decays B^ /}(*)+ should be included 

in the above analysis, since the decay constants Jk and J'k* are also known (the latter can 
be extracted from exclusive r-decay data). On the other hand, Eq. (^) may be used together 
with the experimentally determined value of the QCD coefficient ai to extract yet unknown 
decay constants of spin-one mesons like the ai-meson or the D*-meson without resorting to 
some particular form factor model or to heavy-quark symmetry. 

From our kinematical argument about the equivalence of the lepton pair in the semi- 
leptonic decay and the spin-one particle in the nonleptonic decay it is clear that Eq. (0) 



2We use fg = 205 MeV 
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must be valid, separately, for the different polarizations of the D*-meson in the final state. 
This amounts to the factorization prediction that the polarization of the D*-meson in the 
nonleptonic decay 5° — > D*^ V~ should be equal to the polarization in the corresponding 
semi-leptonic decay at the same q^. This prediction is currently being tested by the CLEO 



collaboration [0]. However, in interpreting the results of such a test, one has to bear in 
mind that in the semi-leptonic as well as in the nonleptonic case the D*-polarization at the 
points = and = q"^^^ is unambiguously determined by kinematics alone to be 100% 
longitudinal and 1/3 longitudinal, respectively. At zero recoil, there is no preferred direction 
and thus the value 1/3 just expresses the fact that there are two transverse, but only one 
longitudinal polarization. At = 0, corresponding to maximum recoil in the semi-leptonic 
decay, the left-handed electron or muon and the right-handed anti-neutrino go off parallel 
to each other, thereby forcing the D* into longitudinal polarization. In the corresponding 
nonleptonic decay S° — D* V~ we know (even without the factorization approximation!) 
that the decay amplitude must be proportional to the polarization vector of the (pseudo- 
) vector meson V. Now, for small = rriy <^ rn^/A the (pseudo-) vector meson V is 
highly relativistic (in the 5-meson rest frame) so that for longitudinal polarization of V 
(and consequently D*) the components of the polarization vector acquire very large values, 
causing longitudinal polarization to dominate. 

The above discussion shows that in comparing polarizations in semi-leptonic and non- 
leptonic decays, one needs polarization data of quite high precision in order to make a 
statement about the validity of factorization. Thus, at low g^, it is the amount of trans- 
verse polarization that has to be measured with a small relative uncertainty. Especially 
for the semi-leptonic decay such a high precision measurement of the g^-dependence of the 
D*-polarization seems hardly possible at present. Fortunately, this seems to be a case where 
heavy-quark symmetry predictions receive only minor corrections. In the infinite quark mass 
limit one finds for the ratio of transverse to longitudinal polarization at some fixed g^ 

dTr ^ 4g^(m| + mg). - g^) 

dTi {mB-mD,Y[{mB + mD'Y - q^] 

which is subject to QCD- as well as l/mg-corrections. The general structure of the correc- 
tions can be found in Ref. fl^. While the Q CD-corrections can be reliably calculated using 



perturbation theory (see e.g. Ref. fl^), the l/mg-corrections are mo del- dependent. We 



have calculated the corrections to Eq. (|3) using the QCD-corrections of Ref. [|18l and the 
l/mg-corrections resulting from a) an analysis of the wave function model of Bauer, Stech 
and Wirbel |1T9| , ^ and b) a sum rule calculation [^. Although corrections to individual 



form factors in both models are as large as 30% at maximum recoil and furthermore vary 
strongly between both models, the correction factor to Eq. (|^) in neither one of the two 
models deviates from one by more than 5% (though the deviations in the two models go 
in opposite directions). In Table |I|, we present the prediction for the D*-polarization in 
semi-leptonic -B-decays as a function of g^ obtained from Eq. ([7p. The quoted errors result 
from the conservative estimate of a 10% relative uncertainty for Tt/'^l at maximum recoil 
(i.e. at g^ = 0), decreasing linearly to the point of zero recoil, where the polarization is fixed 
model-independently. 

In nonleptonic S-decays, the only polarization measurement presently available is that 
of the D^-polarization in the decay 5° D*+ g- . CLEO finds VT/Ttot. = (10 ± 9)% P 
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which has to be compared with the 12% transverse polarization predicted for the semi- 
leptonic decay at = (see Table In order for this test to be sensitive to deviations 
from factorization, the experimental uncertainty will have to be reduced. 

The situation may be more favorable in the decay D*^ D* with predicted 48% 

of transverse polarization, hopefully allowing for a measurement with smaller relative un- 
certainties. Also, it will be particularly interesting to see whether this decay obeys the 
factorization prediction, as this would indicate that the factorization assumption may be 
justified even in decays with only medium energy release. However, one should keep in mind 
that the QCD coefficient Oi drops out of the ratio Tt/T^, so that from polarization tests 
alone it will not be possible to decide whether the short range QCD corrections represented 
by the values of ai and a2 are really independent of the energy release. To this end, one 
would like to test the validity of Eq. with V = D* using a value for the decay con- 
stant of the D* as determined independently from a measurement of the rate for the decay 
Dg ^ liu (employing /^j^ ~ predicted by heavy-quark symmetry). In the absence of 
such a measurement fo* may be taken from sum rule or lattice calculations. 
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TABLES 



TABLE L Determination of the QCD coefficient ai from several nonleptonic B decay modes 
as a test of the factorization assumption. The data are taken from CLEO and ARGUS. The 
theoretical predictions for the branching ratios in the last two rows are those of Ref. @ 



Quantity 


Experiment 


Theory 


Ol 


K [see Eq. (§] 


1.29± 0.22 


0.97af 


1.15±0.10 


R; [see Eq. (|)] 


3.0± 0.7 


2.37af 


1.13±0.13 


p* / n* 


2.5 ±0.6 


f,/f^ = 2.4 






0.28 ± 0.05 


0.214a? 


1.15 ±0.10 


BR(i?o ^ D+g ) 


0.74 ± 0.22 


0.502a? 


1.21 ±0.18 



TABLE IL Amount of transverse polarization (in %) of the D* in semi-leptonic i?-decay. 








2 




ml* 


'Zmax 


drr/dEtot. 





12± 1 


26± 2 


48± 1 


2 
3 
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